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We present results on the mass of the nucleon and the A using two dynamical degenerate twisted 
mass quarks and the tree-level Symanzik improved gauge action. The evaluation is performed at 
four quark masses corresponding to a pion mass in the range of about 300-600 MeV on lattices 
of 2.1-2.7 fm at three lattice spacings less than 0.1 fm. We check for cut-off effects by evaluating 
these baryon masses on lattices of spatial size 2.1 fm at j3 = 3.9 and j3 = 4.05 and on a lattice of 
2.4 fm at f3 = 3.8. The values we find are compatible within our statistical errors. Lattice results 
are extrapolated to the physical limit using continuum chiral perturbation theory. Performing a 
combined fit to our lattice data at /3 — 3.9 and /3 = 4.05 we find a nucleon mass of 963 ± 12(stat.) ± 
8(syst.) MeV where we used the lattice spacings determined from the pion decay constant to convert 
to physical units. The systematic error due to the chiral extrapolation is estimated by comparing 
results obtained at 0{p i ) and 0(p 4 ) heavy baryon chiral perturbation theory. The nucleon mass at 
the physical point provides an independent determination of the lattice spacing. Using heavy baryon 
chiral perturbation theory at 0(p 3 ) we find ap^.g = 0.0889 ± 0.0012(stat.) ± 0.0014(syst.) fm, and 
= 0.0691 ± 0.0010(stat.) ± 0.0010(syst.) fm, in good agreement with the values determined 
from the pion decay constant. Using results from our two smaller lattices spacings at constant 
rom^ we estimate the continuum limit and check consistency with results from the coarser lattice. 
Results at the continuum limit are chirally extrapolated to the physical point. Isospin violating 



lattice artifacts in the A-system are found to be compatible with zero for the values of the lattice 
spacings used in this work. Performing a combined fit to our lattice data at (3 = 3.9 and f3 = 4.05 
we find for the masses of the A++'~ and A+'° 1315 ± 24(stat.) MeV and 1329 ± 30(stat.) MeV 
respectively. We confirm that in the continuum limit they are also degenerate. 



PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 
Keywords: Nucleon mass, A mass, Lattice QCD, chiral effective theories 
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I. INTRODUCTION 

Twisted mass fermions provide an attractive formulation of lattice QCD that allows for automatic 0(a) improve- 
ment, infrared regularization of small eigenvalues and fast dynamical simulations [1, 2, 3, 4, 5]. A particularly 
attractive feature is that automatic 0(a) improvement is obtained by tuning only one parameter requiring no further 
improvements on the operator level. A tree-level analysis of cut-off effects for twisted mass fermions has been pre- 
sented in Ref. [6], while a preliminary non-perturbative investigation on scaling of several observables is carried out 
in Ref. [7] . Recent simulations with two degenerate flavors of dynamical Wilson twisted mass fermions demonstrate 

that pion masses of ~ 300 MeV can be reached using Hybrid Monte Carlo methods [5, 8, 9]. The theoretical 
framework to include the strange and charm quarks has been layed out and practical simulations are being investi- 
gated [10, 11, 12]. Important physical results are emerging using gauge configurations generated with two degenerate 
twisted quarks: In the meson sector very precise results on the pion mass and decay constant led to the determination 
of the low energy constants l 3 , f 4 , F and B [7, 13, 14] to an accuracy that had an immediate impact on chiral 
perturbation theory (%PT) predictions [15]. Accurate results on the pion form factor are obtained [16] using the 
"one-end-trick" method developed in Refs. [17, 18]. The kaon system is studied in a partially quenched approach by 
implementing non-degenerate valence twisted mass quarks maintaining automatic 0(a) improvement [19, 20, 21, 22]. 
After determining the average up and down quark mass and the strange quark mass, the kaon decay constant is 
extracted [23, 24]. In a similar approach first results on the charm quark mass and decay constant are obtained [25]. 
Preliminary results on the first moment of the pion quark distribution function were reported in Ref. [26] . 

In this work we present a detailed analysis of results in the light baryon sector, a subset of which was given in 
Ref. [27]. Using two dynamical degenerate twisted mass quarks we evaluate the mass of the nucleon and A for 
pion masses down to about 300 MeV. We use the tree-level Symanzik improved gauge action [28]. We perform the 
calculation using three different lattice spacings corresponding to (3 = 3.8, [3 — 3.9 and (3 — 4.05 to check cut-off 
effects, where (3 = 6/g 2 with g being the bare coupling constant. For each value of (3 we have configurations at four 
different values of the bare quark mass chosen so that the pion masses are in the range of about 300 MeV to 600 MeV. 
These gauge configurations belong to the same ensembles as those analyzed for the evaluation of the pion mass and 
decay constant. The values of the lattice spacing extracted from the pion decay constant are 0,3=3.8 — 0.0995(7) fm, 
0/3=3.9 — 0.0855(5) fm and 0/3=4.05 = 0.0667(5) fm [13, 14] and will be used in this work. At /? = 3.9, for the smallest 
pion mass, there are gauge configurations at two different volumes enabling us to assess finite volume effects. 

Chiral perturbation theory has been successfully applied in the extrapolation of lattice data obtained with twisted 
mass fermions in the pion sector yielding an accurate determination of the relevant low energy constants. Applying 
xPT to the baryon sector is more involved and several variants exist. However, to leading one-loop order, the result 
is well established and the quality of our lattice results allows for extrapolation to the physical point using this lowest 
order result. Performing a combined fit to our lattice data at (3 = 3.9 and (3 = 4.05 using the leading one-loop order 
result we find a nucleon mass of 963 ± \2{stat.) MeV, where we convert to physical units using the lattice spacing 
determined from / x . We would like to point out that in most other chiral extrapolations of lattice data the physical 
point is included in the fits and therefore such a consistency check cannot be made. The nucleon mass at the physical 
point provides an independent determination of the lattice spacing. We find that the lattice spacing thus determined 
is in good agreement with the value extracted in the pion sector. This is a non-trivial check of our lattice formulation 
and of the smallness of the systematic errors involved. To assess systematic errors due to the chiral extrapolation we 
perform chiral fits to the nucleon and A mass using higher order chiral perturbation theory results, which also include 
explicitly the A degree of freedom. 

One of the drawbacks of twisted mass fermions is the 0{a 2 ) breaking of isospin symmetry, which is only restored 
in the continuum limit. In the baryon sector we can study isospin breaking by evaluating the mass difference between 
A ++ (A~) and A + (A°). Unlike in the pion sector, where disconnected contributions enter in the evaluation of the 
mass of the n°, here there are none. We can therefore obtain an accurate evaluation of isospin splitting and its 
dependence on the lattice spacing. We find no isospin splitting within our statistical accuracy. This is in agreement 
with a theoretical analysis [29, 30] that shows potentially large 0{a 2 ) flavor breaking effects to appear in the 7r°-mass 
but to be suppressed in other quantities. Like in the nucleon case, we perform a combined fit to our lattice data at 
(3 = 3.9 and (3 = 4.05 for the mass of the A ++ '~ and A +, ° using the lowest one-loop order chiral perturbation result. 
We find for the mass of the A++-- and A+'° 1315±24(sioi.) MeV and 1329± 30(stat.) MeV respectively. We confirm 
that in the continuum limit they are also degenerate. 

This paper is organized as follows: In Section II we present our lattice action and in Section III we explain our 
lattice techniques. In Section IV we discuss lattice artifacts and in Section V we give results on the nucleon and A 
mass and also describe the chiral extrapolations. Finally in Section VI we provide a summary and conclusions. 
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II. LATTICE ACTION 



For the gauge fields we use the tree-level Symanzik improved gauge action [28] , which includes besides the plaquette 
term U^ 1 ^ also rectangular (1 x 2) Wilson loops 

S 9 = ^J2\ b " E {l"Re'&(^)}+6 1 J2 {l-ReTr(C/^)}J (1) 

with 6i = —1/12 and the (proper) normalization condition b n = 1 — 8b\. Note that at b\ — this action becomes the 
usual Wilson plaquette gauge action. 

The fermionic action for two degenerate flavors of quarks in twisted mass QCD is given by 

S f = ft4 E ( D w[U] + m + i^ 5 T 3 ) X (x) (2) 

x 

with r 3 the Pauli matrix acting in the isospin space, /i the bare twisted mass and the massless Wilson-Dirac operator 
given by 

I OT 

Dw[u} = - lfi (v^ + v;)--v^; . (3) 



where 

1 



V^(a:) 



a 



UUx)i/>(x + ap,) - i>(x) 



and V*ip(x) 



a 



Ufj,{x — afi)ip(x — afi) — ip(x) 



(4) 



Maximally twisted Wilson quarks are obtained by setting the untwisted quark mass too to its critical value m cr , while 
the twisted quark mass parameter \i is kept non-vanishing in order to work away from the chiral limit. In Eq. (2) the 
quark fields \ are m the so-called " twisted basis" . The " physical basis" is obtained for maximal twist by the simple 
transformation 

tp{x) = cxp ( ^75T 3 j x(x), ip(x) = x(x) exp ( ^75T 3 J • (5) 



In terms of the physical fields the action is given by 



Sf = a 4 J2 fa) (^7,[V, + V;] - z 75 r 3 (-y V,V; + to ct ) + M ) Mx) 



(6) 



In this paper, unless otherwise stated, the quark fields will be understood as "physical fields", ip, in particular when 
we define the baryonic interpolating fields. 

A crucial advantage of the twisted mass formulation is the fact that, by tuning the bare untwisted quark mass Too to 
its critical value to ct , all physical observables are automatically 0{a) improved. In practice, we implement maximal 
twist of Wilson quarks by tuning to zero the bare untwisted current quark mass, commonly called PCAC mass, 
™PCAC, which is proportional to m — m cr up to 0{a) corrections. As detailed in Ref. [31], topcac is conveniently 
evaluated through 

r Ex(d 4 ^(x,t)P b (0)) 
topcac = hm — = = , 0=1,2 , 7 

where A 1 ^ = Xlnls^X IS the axial vector current and P b = Xl'a\X t ne pseudoscalar density in the twisted basis. The 
large t/a limit is required in order to isolate the contribution of the lowest- lying charged pseudoscalar meson state 
in the correlators of Eq. (7). This way of determining topcac is equivalent to imposing on the lattice the validity of 
the axial Ward identity d^A^ = 2topcac-P 6 , b = 1,2, between the vacuum and the charged zero three-momentum 
one-pion state. When too is taken such that topcac vanishes, this Ward identity expresses isospin conservation, as 
it becomes clear by rewriting it in the physical quark basis. The value of m CT is determined at each (5 value at the 
lowest twisted mass, a procedure that preserves 0(a) improvement and keeps 0{a 2 ) small [31, 32]. 

The twisted mass fermionic action breaks parity and isospin at finite lattice spacing, as it is apparent from the 
form of the Wilson term in Eq. (6). In particular, the isospin breaking in physical observables is a cut-off effect 
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P = 4.05, a = 0.0667(5) fm 



32^ x 64, L s = 2.1 fm 



an 0.0030 0.0060 0.0080 0.0120 

m n (GcV) 0.3070(18) 0.4236(18) 0.4884(15) 0.5981(18) 



24 a x 48, L s = 2.1 fm 
32 3 x 64, L s = 2.7 fm 



(3 = 3.9, a = 0.0855(5) fm 



an 0.0040 0.0064 0.0085 0.010 
(GeV) 0.3131(16) 0.3903(9) 0.4470(12) 0.4839(12) 
an 0.0040 
m n (GeV) 0.3082(55) 



24 a x 48, L s = 2.4 fm 



13 = 3.8, a = 0.0995(7) fm 



an 0.0060 0.0080 0.0110 0.0165 

m n (GcV) 0.3667(17) 0.4128(16) 0.4799(9) 0.5855(10) 



TABLE I: The parameters of our calculation. 



of 0(a 2 ) [2]. However the up- and down-propagators satisfy G u (x,y) = ^sGj(y, x)^s, which means that two-point 
correlators are equal with their hermitian conjugate with u- and d-quarks interchanged. Using in addition that the 
masses are computed from real correlators, it leads to the following pairs being degenerate: tt + and 7r~, proton and 
neutron and A ++ (A + ) and A - (A ). A theoretical analysis [29] shows that potentially large 0(a 2 ) effects that appear 
in the 7r°-mass are suppressed in other quantities. Calculation of the mass of 7r°, which requires the evaluation of 
disconnected diagrams, has been carried out confirming large 0(a 2 )-effects. In the baryon sector we study isospin 
breaking by evaluating the mass difference between A ++ (A~) and A + (A°). Since no disconnected contributions 
enter, the mass splitting can be evaluated using fixed source propagators. An accurate evaluation of the isospin 
splitting and its dependence on the lattice spacing is an important component of this work. Examining the size 
of isospin breaking is a crucial aspect in particular regarding future applications of twisted mass fermions to study 
baryon structure. We find that the isospin breaking for the values of the lattice spacing considered in this work is 
consistent with zero within our statistical accuracy. Taking the continuum limit of our lattice results we confirm that 
A ++, ~ and A +, ° are indeed degenerate leading to the same mass at the physical point. 



III. LATTICE TECHNIQUES 

The simulation parameters were chosen such that the pion mass ranges from about 300-600 MeV. The lattice 
volumes and masses used in this calculation are collected in Table I. Finite size effects are examined using the smallest 
pion mass at (3 = 3.9 as finite volume effects are largest. At this mass we have simulations on lattices of spatial size, 
L s ~ 2.1 fm and L s ~ 2.7 fm. 

In order to estimate finite o-effects and the continuum limit we use two sets of results at (3 = 3.9 and (3 = 4.05. 
Although a further set of gauge configurations at (3 = 3.8 is analyzed this set is not used to extrapolate to the continuum 
limit. The reason is that the performance of the HMC algorithm that we use for the simulations deteriorates when 
we go to small n values on this coarser lattice. The long autocorrelation times of the plaquette and the PCAC mass 
that we observe [13] make the tuning to maximal twist less reliable than for the finer lattices. An analysis aimed at 
quantifying the impact of possible numerical errors from the tuning procedure on observables [7] is still in progress. 
Therefore the set at (3 = 3.8 is used only as a cross-check and to estimate cut-off errors. 



A. Interpolating fields 



The masses of the nucleon and the A's are extracted from two-point correlators using the standard interpolating 
fields, which for the proton, the A ++ and A + , are given by 



Jp 


= £abc(ulCj 5 d b )u c , 


tA» IT 
J A++ ~~ e abc\U a 


CY u b)u c , 


J A + 


1 

" vf abc 


2(u T a C^d b )u c 


+ (u T a cru b )d c 


1 



where C = 7472- 

Local interpolating fields are not optimal for suppressing excited state contributions. We apply Gaussian smearing to 
each quark field, g(x, t) [33, 34]. The smeared quark field is given by g smoar (x, t) = J2 y F{ x , Y! U (t))q(y, t) using the 
gauge invariant smearing function 



F(x,y;U(t)) = (l + aH) n ( X ,y;U(t)), 



(9) 



■5 



constructed from the hopping matrix understood as a matrix in coordinate, color and spin space, 

3 



ff(x, y; U(t)) = J2 ( U i( x > t) 5 *,y-ai + U H X - °>h *)<W+a?) . 

<=i ^ ' 



(10) 



The parameters a and n are varied so that the root mean square (r.m.s.) radius obtained using the proton interpolating 
field is in the range of 0.3-0.4 fm. In Fig. 1 we show lines of constant r.m.s radius as we vary a and n. The larger 
the n the more time consuming is the smearing procedure. On the other hand, for a ~ 1, increasing further a does 
not reduce n significantly. Therefore, we choose a value of a large enough so that the weak a-dependence sets in, and 
we adjust n to obtain the required value of the r.m.s radius. We consider two sets for these parameters giving r.m.s 
radius 0.32 fm and 0.41 fm, as shown in Fig. 1. For each set of parameters we evaluate the effective mass as 



p 



\og(Cp(t)/C P (t-l)) 



where Cp{t) is the zero- momentum two-point correlator of the particle P given by 

Cp{t) — Tr(l zb ^4) ^ ^ (</p(x s j n k, t s ink)^P (^sourco ^source)) i ^ — ^sink ~~ ^-source 



(ii) 



(12) 



In Fig. 2, we show the nucleon effective mass, m^g in lattice units for 10 configurations at (3 = 3.9 and a/i — 0.0085. 
For the optimization of the parameters we apply Gaussian smearing at the sink, whereas for the source we use local 
interpolating fields so that no additional inversions are needed when we change a and n. As can be seen, for both 
sets of smearing parameters, the excited state contributions are suppressed with the set a — 4, n — 50 producing 
a plateau a couple of time slices earlier. If, in addition, we apply APE smearing [35] to the spatial links that enter 
the hopping matrix in the smearing function, then gauge noise is reduced resulting in a better identification of the 
plateau. Therefore for all computations at p = 3.9 we use Gaussian smearing with a = 4 and n — 50. Having chosen 
the smearing parameters, for the results that follow, we apply smearing at the source and compute the mass using 
both local (LS) and smeared sink (SS). For (3 — 4.05 we readjust the parameters so that the nucleon r.m.s radius is 
still about 0.4 fm, obtaining a = 4 and n = 70. For (3 — 3.8 we use a = 4 and 30 to keep the r.m.s. radius at the 
same value. Also for these two values of (3 we apply APE smearing to the gauge links that are used in F(x, y; U(t)). 

There are other methods to enhance ground state dominance besides Gaussian smearing. Smearing based on link 
fuzzing has been effectively used in the pion sector. In this work, having optimized our parameters for Gaussian 
smearing we use only local and Gaussian-smeared interpolating fields. 
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FIG. 1: Lines of constant r.m.s radius as a function of 
the smearing parameters a and n. The asterisk shows the 
values a = 2.9, n — 30 and the cross a = 4.0, n = 50. 



FIG. 2: m cff versus time separation both in lattice units. 
Crosses show results using local sink and source (LL), 
circles (asterisks) using Gaussian smearing at the sink 
(SL) with a = 2.9 and n = 30 ( a = 4 and n = 50), 
and filled triangles with a = 4 and n = 50 and APE 
smearing. The dashed line is the plateau value extracted 
by fitting results when APE smearing is used. 
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B. Two point correlators 

The lowest energy state with which the nucleon interpolating field given in Eq. (8) has a non-vanishing overlap is 
the proton state \p(p, s)) 

(0\J p \p(p,s)) = Z p u(p,s). (13) 
Z p is a constant overlap factor and u(p, s), with s € { — 1/2, +1/2}, is a solution of the Dirac equation 

( J /-m N )u = 0. (14) 
Averaging over the nucleon spins and choosing the nucleon rest frame, we are led to the two point correlator 

Cjy(t) = — Tr(l i ^4) ^ ^ (t//v(x s inkj ismk)t/7v(x SO urcc; ^source)); t — tsink ^source ■ (15) 



Space-time reflection symmetries of the action and the antiperiodic boundary conditions in the temporal direction 
for the quark fields imply, for zero three-momentum c 
extracted from the exponential decay of the correlator 



for the quark fields imply, for zero three-momentum correlators, that C^(t) = —C N (T — t). The nucleon mass is 



C N (t) = C+(t)-C„(T-t). (16) 
To increase the precision we also average over the proton and neutron correlators since these are degenerate in mass. 
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n Without projection 
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FIG. 3: Comparison of effective masses for A ++ ' for 
aii = 0.0085 at (3 = 3.9 on the lattice volume 24 3 x 48, 
obtained with (filled triangles) or without (open squares 
shifted to the left for clarity) spin projection, using a 
sample of 90 configurations. The mass difference with 
projection and without projection is much smaller than 
the statistical error. 
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FIG. 4: Comparison of effective masses for A + for 
a/i = 0.011 at P = 3.8 on the lattice volume 24 3 x 48, 
obtained with 3/2-spin (filled triangles) or with 1/2-spin 
projection, using a sample of 50 configurations. 



In analogy to Eq. (13), the overlap of the A + interpolating field with the A + state is given by 

(0|J£ + |A+( P , S )) =Z A+ u»(p,s). 
Every vector component of the Rarita-Schwinger spinor satisfies the Dirac equation 

(y-m A K = // = 1...4, 

and in addition the auxiliary conditions 

= and = 



(17) 
(18) 
(19) 



are fulfilled. The four independent solutions are labeled by s S {—3/2,-1/2,1/2,3/2}. The A interpolating fields 
as defined in Eq. (8) have overlap also with the heavier spin-1/2 excitations. These overlaps can be removed when 



7 



the conditions in Eq. (19) are enforced on the interpolating fields. This can be achieved by the incorporation of a 
spin-3/2 projector in the definitions of the interpolating fields 

^3/2 = ^3/2 J" A ( 2 0) 

= ^ - 1^1" - &W + V^f V)- (21) 

Similarly the spin 1/2-interpolating field, Jy 2 , that has only overlap with the 1/2-state, is obtained by acting with the 
spin 1/2-projector P^J 2 = 9 tiu — P3J2 011 ^a- Using any of the three interpolating fields, the A masses are extracted 
from the two-point functions 

1 3 

^a(^) gTr(l it ^4) ^ ^ ^ ^ (^A (-X-sinki tsinkjJ A ( x sourcc ; ^source)) ? t ^sink ^source- (22) 

x aink i=l 

Fig. 3 compares effective masses extracted from correlation functions with and without the spin 3/2 projection at 
(3 = 3.9. For this comparison we use 90 configurations, a number sufficient for the purpose of this check. The results 
for the effective mass are hardly affected by including the spin-3/2 projector even at very short time separations. 
This is because the overlap of the interpolating field with the spin- 1/2 state is small, a property that holds at all 
values of (3 This is clearly seen in Fig. 4 at (3 = 3.8 where the effective mass obtained using the spin-1/2 projected 
interpolating field Jy 2 is much more noisy than with Jy 2 due to the small overlap with the spin-1/2 state. This 
behavior is in agreement with the results of Rcf. [36] where the same spin projections were implemented. Since 
the impact on the plateau value is negligible compared to the statistical uncertainty, we use only the non-projected 
interpolating fields from here on. We average the correlators of A ++ with A~ as well as A + with A . We do not 
average the A ++ and A + correlators as they differ by an 0(a 2 ) isospin breaking effect. 



C. Effective masses 
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24 3 


x 48 








0.0040 


LL 


471 




0.13587(68) 


0.511(11) 


0.699(15) 


0.708(25) 


0.0040 


LS 


419 




0.13587(68) 


0.521(6) 


0.694(11) 


0.717(13) 


0.0040 


ss 


419 




0.13587(68) 


0.515(5) 


0.682(12) 


0.697(17) 


0.0064 


LL 


199 




0.16937(36) 


0.565(11) 


0.727(15) 


0.763(15) 


0.0064 


LS 


235 




0.16937(36) 


0.565(6) 


0.715(13) 


0.742(7) 


0.0064 


SS 


235 




0.16937(36) 


0.561(4) 


0.710(11) 


0.711(10) 


0.0085 


LL 


153 




0.19403(50) 


0.568(8) 


0.754(12) 


0.776(22) 


0.0085 


LS 


186 




0.19403(50) 


0.581(6) 


0.746(9) 


0.751(12) 


0.0085 


SS 


186 




0.19403(50) 


0.580(6) 


0.738(11) 


0.742(9) 


0.0100 


LL 


173 




0.21004(52) 


0.613(6) 


0.823(7)* 


0.767(12) 


0.0100 


LS 


213 




0.21004(52) 


0.595(7) 


0.742(7) 


0.760(7) 


0.0100 


SS 


213 




0.21004(52) 


0.589(9) 


0.750(10) 


0.755(10) 








32 a 


x 64 








0.0040 


LS 


201 




0.13377(24) 


0.518(5) 


0.672(9) 


0.670(14) 


0.0040 


SS 


201 




0.13377(24) 


0.510(5) 


0.660(9) 


0.660(14) 



TABLE II: Results for the nucleon and A mass at (5 = 3.9 for lattices of size 24 x 48 and 32 x 64. LL stands for local sink 
and local source, LS for local sink and smeared source and SS for smeared sink and smeared source. The results for the pion 
mass are taken from Table 2 of Ref. [13] computed using more gauge configurations than we used for the evaluation of the 
baryon masses as well as a different smearing and therefore are the same for LL, LS and SS. With an asterisk we mark results 
for which the effective mass does not show a good plateau. 

In Fig. 5 we show the nucleon effective masses at /3 = 3.9 on a volume 24 3 x 48 for all the values of n considered. 
We smeared the source as described in the previous section and either use a local sink or smear the sink with the same 
smearing used for the source. As expected, the effective masses are consistent for both smeared and local sink yielding 
asymptotically the same constant. We fit the effective mass to a constant in the region where m e s(t) becomes time 
independent (plateau region) and vary the lower t-range of the fit so that x 2 per degree of freedom (d.o.f.) becomes 
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FIG. 5: Nucleon effective mass (LS: asterisks, SS: open 
triangles) for — 3.9 versus time separation in lattice 
units, for a/t = 0.010 (upper left), 0.0085 (upper right), 
0.0064 (lower left) and 0.0040 (lower right). The constant 
lines are the best fits to the data over the range spanned 
by the lines. 
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FIG. 6: A ++ '~ (asterisks) and A + '° (open triangles) ef- 
fective masses using smeared-smeared (SS) correlators 
for f3 = 3.9 versus time separation in lattice units, for 
an = 0.010 (upper left), 0.0085 (upper right), 0.0064 
(lower left) and 0.0040 (lower right). The straight lines 
are the best fits to the data over the range spanned by 
the lines. The solid line is for A ++ '~ and the dashed for 
A + '° and they coincide. 



less than one. We take this value for the mass of the nucleon. In Fig. 6 we show, for the same //-values, effective 
masses for the A ++, ~ and A + '° using smeared source and sink. We fit in the same way as in the nucleon case to 
extract the mass of the A. As can be seen, the quality of the plateaus in the nucleon case is better than in the case 
of the A. This explains why results on the A mass have larger errors. The errors are evaluated using jackknife and 
the T-method [38] to check consistency. The integrated auto-correlation times for our baryonic observables are very 
small for our configuration ensembles. Since for our computation we use gauge configurations that are separated by 
8-20 trajectory lengths, autocorrelations are negligible. 

The resulting masses using local and smeared interpolating fields are summarized in Tabic II for (3 — 3.9 using 
lattice sizes of 24 3 x 48 and 32 3 x 64, while those obtained for (3 — 4.05 on a lattice volume of 32 3 x 64 are reported 
in Table III. Results obtained at j3 = 3.8 are given in Table IV. The mass of the pion listed in Table II is taken from 
Ref. [13] and was evaluated using a larger set of configurations applying a different smearing than the one used in this 
work. A detailed description of this evaluation as well as the error analysis is presented in Ref. [31]. The pion masses 
given in Tables III and IV are again obtained in a separate evaluation [37]. 



IV. LATTICE ARTIFACTS 



A. Finite volume effects 



At (3 — 3.9 and for a// = 0.004 we have gauge configurations on two lattices of different volume. This is the 
smallest //-value considered at (3 — 3.9 and it is the one that potentially can have the largest finite volume effects. 
On the lattice of spatial extension L s — 24 the other three larger //-values satisfy the condition m^Lg > 4, whereas 
for a\i = 0.004 we have m^Ls ~ 3.2. On the 32 3 lattice at afi — 0.004 we have m^Lg > 4. Applying the re-summed 
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SS 
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0.446(6) 
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LL 
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0.1651(5) 


0.465(6) 


0.627(6) 0.637(7) 


0.0080 


LS 
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0.1651(5) 


0.469(4) 


0.590(7) 0.585(9) 


0.0080 


SS 


256 


0.1651(5) 


0.465(5) 


0.594(7) 0.594(8) 
0.670(4)* 0.677(5) 


0.0120 


LL 


157 


0.2025(6) 


0.520(5) 


0.0120 


LS 


162 


0.2025(6) 


0.509(4) 


0.616(7)* 0.623(7) 
0.616(7)* 0.620(7) 


0.0120 


SS 


162 


0.2025(6) 


0.515(3) 



TABLE III: Results for the nucleon and A mass at f3 — 4.05 for the 32 3 x 64 lattice. LL stands for local sink and local source, 
LS for local sink and smeared source and SS for smeared sink and smeared source. The results for the pion mass are computed 
using more gauge configurations than we used for the evaluation of the baryon masses as well as a different smearing [37] and 
therefore are the same for LL, LS and SS. . With an asterisk we mark results for which the effective mass does not show a 
good plateau. Empty entries are due to the absence of a sufficient plateau region. 
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LL 
SL 
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SL 
LL 
SL 
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249 
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0.1852(10) 
0.1852(10) 
0.2085(8) 
0.2085(8) 
0.2424(5) 
0.2424(5) 
0.2957(5) 
0.2957(5) 



0.623(20) 
0.637(9) 
0.676(11) 
0.665(9) 
0.700(9) 
0.699(8) 
0.759(7) 
0.744(8) 



0.792(25) 
0.818(11) 
0.859(11) 
0.827(17) 
0.861(13) 
0.854(14) 
0.948(12) 
0.942(13) 



0.815(28) 
0.824(13) 
0.847(30) 
0.856(24) 
0.893(22) 
0.875(16) 
0.886(25) 
0.946(12) 



TABLE IV: Results for the nucleon and A mass at (3 = 
Table III with SL being a smeared sink and local source. 



3.8 for the 24 3 x 48 lattice. The notation is the same as that of 



Liischer formula to the nucleon mass and using the knowledge of the ttN scattering amplitude to 0{p 2 ) and 0(p 4 ), 
the volume corrections are estimated to be about 3% to 5% [39] for L s ~ 2 fm and m v <~ 300 MeV. In Table II we 
give the results for the nucleon mass using our two lattice volumes. The smaller lattice volume has spatial length very 
close to the 2 fm length of Ref. [39] namely L s ~ 2.1 fm. The results for mjy do not change within our statistical 
accuracy when we use the larger lattice size of L s <~ 2.7 fm. We make the assumption that for the larger lattice 
finite volume corrections have become negligible and take them to be a good approximation to the infinite volume 
results. In other words we take miq{L s = oo) ~ miq{L s = 2.7 fm). This assumption was shown to be valid in the 
pion sector where a finite size analysis was carried out [31]. We define the ratio Rn = AmN/m]y(L s — oo), where 
Am;v = toat(L s = 2.1 fm) — rriN(L s — oo) and estimate Rn with results obtained on our two volumes for the smallest 
pion mass. This gives us an estimate for our finite volume errors. Using the results tabulated in Table II at ay, = 0.004 
we conclude that Amjv is positive as expected. This is also true for the corresponding difference for the masses of 
A ++, ~ and A +, °. In Table V we give the ratios Rn, R a++ - and R&+,o. For the nucleon this ratio is compatible with 
zero and within our accuracy it can be at the most 2%. For the A ++,_ where the statistical errors are smaller than 
for the A + '°, the volume corrections range from 1% to 5%. From this study we conclude that finite volume effects on 
the nucleon mass are negligible whereas for the A we can at most have corrections on the few percent level. 



ay 


Interpolating field 


-Rjv -Ra++.- -Ra+.° 


0.0040 


LS 


0.006(15) 0.033(22) 0.070(26) 


0.0040 


SS 


0.010(14) 0.033(23) 0.056(33) 



TABLE V: Finite volume dependence at /3 = 3.9 for ay — 0.004. Results with a lattice of size 24 3 x 48 are compared to those 
obtained with a lattice size of 32 x 64. For a hadron state P we define Rp = (mp(L a = 2.1 fm) — mp(L s = 2.7 im))/mp(L s = 
2.7 fm) ~ (mp(L s — 2.1 fm) — mp(oo))/mp(oo) assuming the masses at 2.7 fm to be close enough to the infinite volume limit. 
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B. Isospin breaking 
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FIG. 7: The A ++ '" (upper graph) and A+'° (lower 
graph) mass as a function of mjj. for /3 = 3.9 on a lat- 
tice of size 24 3 x 48 (filled triangles) and on a lattice of 
size 32 3 x 64 (open triangles). Results at (3 = 4.05 are 
shown with the filled squares. The physical A mass is 
shown with the asterisk. 
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FIG. 8: The mass splitting between A + '° and A ++ " 
normalized with the mean value of their mass m,A as a 
function of m% in lattice units. Results at (3 = 3.8 are 
shown with the asterisks. The rest of the notation is the 
same as in Fig. 7. 



One of the main goals of this work is to examine isospin breaking in the baryon sector due to lattice artifacts. 
As already explained the proton and the neutron are degenerate. Isospin breaking in the light baryon sector can be 
examined for the A. In Fig. 7 we show results for the mass of A ++ -~ as well as for the mass of A + '°. Results at 
(3 = 3.9 and (3 — 4.05 fall on the same curve pointing to small cut-off effects. Small finite volume effects are visible at 
the smallest pion mass at (3 = 3.9 as discussed in the previous subsection. To check for isospin breaking we plot the 
mass difference between the pairs A ++ ', A~ and A + , A° normalized by the mean value of their mass in Fig. 8 for 
(3 = 3.8, (3 = 3.9 and (3 = 4.05. As can be seen, the splitting is consistent with zero for these values of (3, indicating 
that isospin breaking in the A system is small. 



V. CHIRAL EXTRAPOLATION 



A. Nucleon mass 



We show our results for the nucleon mass as a function of in Fig. 9. The masses are extracted in lattice units. 
To convert to physical units we need to know the value of the lattice spacing. A standard procedure is to match 
the experimental value of to the one obtained on the lattice extrapolated to the physical pion mass. This gives 
a = 0.0855 fm at (3 = 3.9 and a = 0.0667 fm at (3 = 4.05 [14]. We use these values to convert lattice results to physical 
units. The results at these two /3-values fall on a common curve indicating that cut-off effects are small for these 
values of the lattice spacing. In Fig. 9 we include, for comparison, results obtained with dynamical staggered fermions 
from Ref. [40, 41]. Results using these two formulations are consistent with each other. As we already discussed, 
results obtained on lattices of spatial length L s = 2.1 fm and L s = 2.7 fm at (3 = 3.9 for the lowest pion mass are 
consistent indicating that finite volume effects are small for the pion masses used in this work. Therefore as a first 
analysis of our lattice results, we use continuum chiral perturbation theory in an infinite volume to perform the chiral 
extrapolation to the physical point. An analysis carried out after taking the continuum limit will serve as a check of 
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FIG. 9: The nucleon mass as a function of for j3 = 3.9 
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FIG. 10: The variation of the fit parameters m% and ci. 
The elliptical boundary is determined by changing these 
parameters so that the minimal value of \ 2 changes by 
one. The most elongated ellipse is for f3 = 4.05 using 
the nucleon mass at the three lighter pion masses, the 
intermediate is for j3 = 3.9 using all five points and the 
smallest is for a combined fit to both /3-values using a 
total of eight points. 



cut-off effects. The leading one-loop result in heavy baryon chiral perturbation theory (HB^PT) [42] is well known: 

m N = m° N - 4 Cl m2 - -^m^ (23) 

with mP N , the nucleon mass in the chiral limit, and c\ treated as fit parameters. This lowest order result is the 
same in HB^PT with dimensional and infra-red regularization as well as when the A degree of freedom is explicitly 
included. It is also the same in manifestly Lorenz-invariant formulation with infrared regularization. Therefore 
we will use this well established result to predict the nucleon mass at the physical point as well as fix the lattice 
spacing using the experimental nucleon mass as input. The higher order results will only be used to estimate the 
systematic error associated with the chiral extrapolation. We take for f n and 17,4 their physical values, namely 
/ w = 0.092419(7)(25) GeV and gA — 1.2695(29), which is what is customarily done in chiral fits to lattice data on the 
nucleon mass [43, 44, 45]. We will take the experimental values for / w and gA also when using higher order results. 
In higher orders new low energy constants enter, and we also fix their values from experimental data. In order to 
determine the errors on the fit parameters we allow for a variation in the parameters that increases the minimal value 
of x 2 by one. In Fig. 10 we show the boundary of the allowed variation of the parameters. As expected when the 
number of available lattice results increases the error decreases. In Fig. 11 we show fits to the 0(p 3 ) result with the 
error band determined by the maximum allowed variation in the parameters that increase the minimal % 2 by one. 
As can be seen, this 0(p 3 ) result provides a very good fit to our lattice data both at (3 = 3.9 and (3 = 4.05. Since 
finite volume effects are small we use in the fit data on both volumes at (3 = 3.9. In Table VI we give the values 
of the parameters mP N and c\. In this determination we use the lattice spacing determined from f^. In the case of 
(3 = 3.9 we include the result obtained using the larger volume. We would like to stress that, despite the fact that the 
physical point is not included in the fit as customary done in other chiral extrapolations of lattice data. The value of 
the nucleon mass that we find at the physical pion mass using data at both (3 = 3.9 and (3 = 4.05 to fit to the 0(p 3 ) 
HBxPT of Eq. (23) is 963(12) MeV, where the error is only statistical. 
The nucleon sigma term is defined by 

q—u.d ^ 

where we have neglected contributions from other quarks. Following Ref. [46] we use the relation to 2 ~ /1 to evaluate 
on by computing m 2 d J^Ji . Using the value of c\ determined from the nucleon fit we find at the physical point 
(Tat = 66.7 ± 1.3 MeV, where the error is statistical. This value is larger than the prevailing value of 45 ± 8 MeV [47] 
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FIG. 11: On the left set of graphs we show chiral fits to 
the nucleon mass for f3 = 3.9 using 0,3=3.9 = 0.0855 fm 
to convert to physical units. On the right set of graphs 
we show the corresponding chiral fits for f3 — 4.05 using 
a /3=4.05 = 0.0667 fm. The upper most graph shows the 
fit to 0{p i ) HBxPT where we use our results at the three 
lowest values of the pion mass. For the higher order 
fits we perform a simultaneous fit to both f3 = 3.9 and 
[3 — 4.05 always excluding at j3 = 4.05 the result at 
the largest pion mass. The physical point, shown by the 
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but in agreement with a new analysis [48] that includes additional data. Our current calculation does not include a 
dynamical strange quark and a better understanding of this term could come when simulations with dynamical strange 
quarks are available [49] . Note that given the role of the sigma term for what concerns the chiral extrapolation as 
well as its implication in dark matter detection [50] it is clear that a serious effort to better fix its experimental value 
is highly desirable. 

Chiral corrections to the nucleon mass are known to 0(p 4 ) within several expansion schemes. In HBxPT to 0(p 4 ) 
with dimensional regularization [43, 51, 52] the result is given by 
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We take the cut-off scale A = 1 GeV and fix the dimension two low energy constants c 2 = 3.2 GeV -1 [53] and 
C3 = —3.45 GeV -1 [45, 52]. This value is consistent with empirical nucleon-nucleon phase shifts [54, 55]. The 
counter-term E\ is taken as an additional fit parameter. HBxPT with dimensional regularization is in agreement 

o 2 5 

with covariant baryon xPT with infrared regularization up to a recoil term given by g ^™" ^ that is of no numerical 

25D7T / N 

significance [45]. We have included this term in our fits. In the so called small scale expansion (SSE) [45], the A- 
degrees of freedom are explicitly included in covariant baryon x?T by introducing as an additional counting parameter 
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FIG. 13: Continuum limit of the nucleon mass using the 
lowest order HB^PT to interpolate except for the heav- 
iest pion mass where we used linear interpolation. For 
ro and o we use the values determined in the pion sec- 
tor, namely r /a = 4.46(3), (ro = 0.444(4) fm), r /a = 
5.22(2), (r = 0.446(3) fm) and r /a = 6.61(3), (r = 
0.441(4) fm) at (3 = 3.8, (3 = 3.9 and p = 4.05 respec- 
tively [7]. 
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the A-nucleon mass splitting, A = iiia — mN, taking C(A/tojv) ~ Oirii^jm^'). In SSE the nucleon mass is given by 

3(<?i + 3ci) 4 (3^ + 10(2), 
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where R(m„) = \Jm\ 



A 2 



- 1 for 



< A. We 



I 2 cos" 1 (£\ if m v > A and R{m^) = ^A 2 - m 2 log - 

take ca = 1.127 [45], A = 1 GeV and fit the counter-term E\, A different counting scheme, known as 5-scheme, takes 
A/mjy ~ 0(5) and m^/rriN ~ C((5 2 ) [44]. Using the S— scheme in a covariant chiral expansion to order p 4 /A one 
obtains an expansion that has a similar form for the nucleon and A mass. The nucleon mass is given by 
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The ttN and 7rA loop function Vj 00 p is given in Ref. [44] and, following the same reference, we take the value of 
%a — 2.85. Here we use the variant of the (5-scheme that includes the 7rA-loop and adds the fourth order term cim% 
as an estimate of higher order effects, since the complete fourth order result is not available. The parameter c 2 is to be 
determined from the lattice data. The fits using these different formulations are shown in Fig. 11. At j3 = 3.9 shown 
in the left panel we used 0,3=3.9 = 0.0855 fm to convert to physical units. We have four values for the the 24 3 x 48 
lattice and one for the larger lattice. At j3 = 4.05 we only use results at the three smallest pion masses since including 
the result at the largest pion mass yields fits with unacceptably large x 2 /d.o.f.. Therefore only at lowest order x?T 
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C(p 4 /A) 1 5-scheme 


0.886(14) 
0.875(50) 
0.884(51) 
0.867(54) 


-1.21(2) 
-1.23(17) 
-1.19(14) 
-1.29(18) 


0.901(37) 
0.929 
0.944 
0.927 


-1.18(4) 0.889(13) 
-1.10 0.881(42) 
-1.02 0.891(47) 
-1.13 0.873(46 ) 


-1.20(2) 
-1.22(12) 
-1.17(15) 
-1.28(15) 


0.904(14) 
0.893(47) 
0.903(52) 
0.886(51) 


-1.19(2) 0.898(9) 
-1.21(12) 0.889(25) 
-1.16(15) 0.901(30) 
-1.27(16) 0.883(29) 


-1.19(1) 
-1.21(7) 
-1.15(9) 
-1.26(9) 




ml 


Cl 


mi 


ci m u A 


Cl 




ci 


Cl 


A ++,- 


0(p A ) HBxPT 1.248(31) 
0(p 4 /A),5-scheme 1.258(126) 


-1.19(4) 
-1.15(43) 


1.222(68)* 


-1.20(5)* 1.241(27) 
1.347(90) 


-1.21(4) 
-0.85(30) 


1.274(33) 
1.267(80) 


-1.17(4) 1.251(16) 
-1.16(20) 1.261(54) 


-1.20(2) 
-1.16(17) 


A +,u 


0{p A ) HBxPT 1.255(40) 
0(p 4 /A)(5-scheme 1.302(43)* 


-1.20(5) 
-1.03(7)* 


1.261(6)* 


-1.19(8)* 1.256(33) 
1.372(104) 


-1.20(4 ) 
-0.81(32) 


1.264(32) 
1.373(65) 


-1.18(4) 1.262(19) 
-0.85(17) 1.267(42) 


-1.19(3) 
-1.16(12) 



TABLE VI: Fit parameters rn% and in GeV and ci in GeV -1 . Results with an asterisk have x 2 /d.o.f. larger than one. All 
fits to the continuum results excluded the largest value of rom, with the exception of the A + '° in the 5-scheme where to obtain 
a good fit we use all six points. For (3 — 3.9 we use all masses including the results at the larger volume whereas for (3 = 4.05 
we use results at the three smaller pion masses. For the fits to continuum results we give two sets of results: the first set is 
obtained when using linear interpolation to the reference pion masses and the second using C(p 3 ) HBxPT for interpolation. 



(3 = 3.9 (3 = 4.05 P = 3.9 and (3 = 4.05 Continuum 

ra'jy "/3=3.9 m° N a/3=4.05 m% ap=3.9 Q/3=4.05 m% rp 

0{p A ) HBxPT 0.865(2) 0.0886(18) 0.868(5) 0.0708(37) 0.866(1) 0.0889(12) 0.0691(10) 0.868(2) 0.473(9) 
0(p 4 ) HBxPT 0.862(9) 0.0869(46) 0.871 0.0717 0.863(4), 0.0875(26) 0.0681(24) 0.863(8) 0.461(23) 
0(p 4 ) SSE 0.865(10) 0.0876(47) 0.876 0.0724 0.866(9) 0.0884(40) 0.0687(31) 0.866(9) 0.464(24) 

Q(p 4 /A),5-scheme 0.859(11) 0.0865(65) 0.870 0.0717 0.861(9) 0.0873(39) 0.0678(31) 0.861(10) 0.458(25) 



TABLE VII: Determination of the lattice spacing in fm and m% in GeV using the nucleon mass. Fitting to the nucleon 
continuum results obtained by linear interpolation at the five lighter reference pion masses we extract the continuum value of 
ro in fm by constraining the fits to reproduce the physical nucleon mass. 



where we have only two fitting parameters we can perform a fit. For the higher order we give the values of m° N and c\ 
of the curves that pass through all the lattice points. Since cut-off effects are consistent with zero for these two values 
of (3 we can use these two sets of results in a combined fit. For the lattice data at (3 = 4.05 we use 0/3=4.0 = 0.0667 fm 
determined from f n , to convert to physical units. The experimental value of the nucleon is shown with the asterisk. 
In Table VI we give the values of m° N and c\ when simultaneous fits to both j3 = 3.9 and (3 = 4.05 data are done. 
We use the lattice spacings determined from the pion decay constant to convert to physical units [14]. These fits are 
shown in Fig. 11 when using higher order xPT. All formulations provide a good description of the lattice results and 
yield a nucleon mass at the physical point consistent with the experimental value. As already discussed, the value 
of the nucleon mass that we find using Eq. (23) is 963(12) MeV. The corresponding value using 0(p 4 ) HBxPT is 
955(33) MeV. We take the difference between these two mean values as an estimate of the systematic error due to 
the chiral extrapolation and quote 963 ± 12(stat.) ± 8(syst.) MeV as our prediction of the nucleon mass. Within the 
statistical and estimated systematical uncertainty this value is close to the experimental one. Furthermore the values 
extracted for the nucleon at the chiral limit m? N as well as c\ are in agreement in all formulations. In addition the 
value of nucleon ctn term defined in Eq. (24) can be evaluated using HBxPT to 0(p 4 ) of Eq. (25). If we use the next 
to leading order relation between and the quark mass \i q [13, 14, 17] instead of the leading order relation cx [i q 
we find a value of 67 ± 8.0 MeV at the physical point, which is consistent with the value obtained to 0(p 3 ) albeit 
with a larger error. We note that the impact on cat of using Eq. (24) with the next to leading order result, rather 
than the lowest oder relation, between m\ and [i q is small and yields a relative decrease of its value at the physical 
point of about 2% only. 

The consistency between the 0(p 3 ) result and the higher order expansions allows for an extrapolation to the physical 
point and a determination of the lattice spacing using the nucleon mass. Fixing the lattice spacing from the nucleon 
mass allows for a comparison with the value obtained from the pion sector and provides a non-trivial check of our 
lattice formulation. We consider 0,3=3.9 and 0^3=4.05 as independent fit parameters in a combined fit of data at = 3.9 
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r rn n 


roiriN 


rom A ++,- 


rorn A +,o 


P = 3.8 



0.70 
0.80 
0.90 
1.00 
1.10 
1.25 



0.70 
0.80 
0.90 
1.00 
1.10 
1.25 



0.70 
0.80 
0.90 
1.00 
1.10 
1.25 



0.70 
0.80 
0.90 
1.00 
1.10 
1.25 



2.654(95 
2.804(45 
2.935(40 
3.044(45 
3.133(34 
3.256(28 



2.666(27 
2.840(27 
2.947(22 
3.013(31 
3.132(24 
3.313(76 



2.676(56 
2.843(71 
2.903(40 
3.003(40 
3.086(24 
3.287(24 



2.667(24 
2.840(25 
2.936(19 
3.009(25 
3.109(17 
3.289(23 



2.668(55 
2.807(41 
2.933(27 
3.043(21 
3.129(31 
3.201(63 



2.672(17 
2.809(13 
2.933(12 
3.041(14 
3.125(20 
3.194(33 



2.665(47 
2.795(36 
2.910(24 
3.007(16 
3.079(19 
3.126(42 



2.671(16 
2.807(12 
2.929(11 
3.025(11 
3.101(14 
3.168(26 



3.596(119) 
3.637(54) 
3.667(75) 
3.731(83) 
3.831(63) 
4.086(62) 



3.511(69 
3.614(47 
3.701(28 
3.767(31 
3.804(57 
3.793(111) 



P = 3.9 
3.548(60) 3.481(35 



3.621(65) 
3.720(56) 
3.841(60) 
3.917(50) 
4.028(66) 



3.614(25 
3.734(16 
3.837(18 
3.916(30 
3.976(56 



P = 4.05 



3.565(124) 
3.627(137) 
3.669(45) 
3.771(45) 
3.932(43) 
4.017(48) 



3.428(81 
3.571(59 
3.701(35 
3.815(18 
3.905(35 
3.982(85 



Continuum 



3.551(54) 
3.622(59) 
3.687(35) 
3.794(36) 
3.926(33) 
4.021(39) 



3.472(32 
3.608(23 
3.728(15 
3.826(13 
3.916(23 
3.978(47 



3.502(161 
3.641(65 
3.790(111 
3.885(120 
3.918(72 
4.127(64 



3.632(88 
3.641(61 
3.725(52 
3.861(46 
3.943(50 
4.061(105 



3.550(97 
3.627(107) 
3.690(45 
3.793(46 
3.936(50 
4.038(49 



3.595(65 
3.638(53 
3.705(34 
3.827(33 
3.939(35 
4.042(44 



3.528(79) 
3.652(53) 
3.760(30) 
3.850(33) 
3.915(64) 
3.950(127) 



3.498(46) 
3.633(32) 
3.754(20) 
3.859(17) 
3.938(31) 
4.000(62) 



3.477(71) 
3.609(51) 
3.727(30) 
3.827(18) 
3.902(37) 
3.953(84) 



3.492(38) 
3.626(27) 
3.746(16) 
3.844(12) 
3.924(24) 
3.984(50) 



TABLE VIII: Results for the nucleon, A ++ '~ and A + '° mass interpolated at the same value of rom^ for the three P values. 
The continuum limit is then taken at constant rom^ using the results at p — 3.9 and p — 4.05. We give results using linear 
interpolation in second, fourth and sixth columns whereas in the third, fifth and seventh columns we give the results using 
lowest order HBxPT. 



and P = 4.05 where the physical nucleon mass is included with no error. In this way the lattice spacings can be 
determined solely by using as input the nucleon mass at the physical point. Using the leading one- loop result we find 
0^=3.9 = 0.0889(12) fm and (1/3=4.05 = 0.0691(10) fm. The quality of these fits are shown in Fig. 12. The values of 
the lattice spacing obtained to 0(p 4 ) using Eq. (25) are given in Table VII. Both SSE and the S— scheme defined by 
Eqs. (26) and (27), which include explicitly A-degrees of freedom, yield values that are consistent with those obtained 
in HBxPT. The variation in the value of a in the different chiral extrapolation schemes gives an estimate of the 
systematic error involved in the chiral extrapolation. We take the difference between the mean values obtained using 
0(p 3 ) and 0(p 4 ) HBvPT as an estimate of the systematic error. Our lattice spacings fixed using the nucleon mass are 
therefore ap= 3 . 9 = 0.0889 ± 0.0012(stat.) ± 0.0014(sysi.) fm and 0^=4.05 = 0.0691 ± 0.0010(siai.) ± 0.0010(sysi.) fm 
and are on the upper bound of the values obtained using f n . 

The physical spatial volumes of the 24 3 lattice at /3 = 3.9 and that of the 32 3 lattice at P = 4.05 are about 
(2.1) 3 fm 3 . Bearing in mind that volume corrections for this lattice size are shown to be small we use results obtained 
on these two almost equal volumes to estimate our masses at the continuum limit. In order to take the continuum 
limit we interpolate data, expressed in units of ro, at the same value of r^m^, where the Sommer parameter ro is 
determined from the force between two static quarks. We use ro/a = 4.46(3), r /a = 5.22(2) and r /a = 6.61(3) for 
P = 3.8, P = 3.9 and /3 = 4.05 respectively [7]. The values of r^m^ that we choose are close to the pion mass values 
where our computation is done and are given in Table VIII. We use a linear interpolation or the fit curves determined 
using chiral effective theories to obtain the nucleon mass at these reference values of r^m^. This procedure is done 
for our three p- values. We use the results at constant r m^ at p = 3.9 and P = 4.05 to estimate the continuum limit 
by fitting to a constant under the assumption that residual cut-off effects on the pion and nucleon masses as well as 
on ro at p = 3.9 and = 4.05 are negligible. This assumption is corroborated by our lattice data shown in Fig. 9 
Results at P = 3.8 at the same value of r^m^ serve as a check for the consistency of this procedure. This is illustrated 
in Fig. 13 where results at P = 3.8 are consistent with the constant fit. The results at the continuum limit are then 
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chirally extrapolated. The parameters obtained are given in Table VI and the fits are shown in Fig. 14 where we 
excluded the heaviest pion mass from these fits. The values of wP N and c\ obtained from the fits to continuum results 
are consistent with the values obtained using results at finite a. This demonstrates that cut-off effects are small. 

The value of the parameter ro can be determined from our results in the continuum limit using the value of the 
physical nucleon mass. We give the extracted values in Table VII where we used linear interpolation to obtain 
the nucleon mass at the reference values of r^m,^. The values extracted for ro in the continuum limit using these 
fits are consistent. Had we used chiral fits at 0(p 3 ) to interpolate the value extracted would change by 0.004 fm 
and at 0{p 4 ) by 0.002 fm. These chan ges are smaller than the statistical errors. We again take the variation in 
the value of ro at 0(p 3 ) and 0(p 4 ) HB^PT as an estimate of the systematic error due to the chiral extrapolation. 
Using the values given in Table VII this difference is 0.012 fm. We add to this error the variation in the values 
obtained using a linear interpolation scheme and the 0(p 3 ) fit, which is 0.004 fm. Therefore the value that we find 
is r = 0.473 ± 0.09(stat.) ± 0.016(syst.) fm. This value of r , like for the lattice spacing, is at the upper bound of 
the value ro = 0.444(3) fm [14] determined using f v . The validation of these consistency checks suggests that lattice 
artifacts that can affect the value of the lattice spacing when using different observables are small. 



B. A mass 



We perform a similar analysis as for the nucleon mass in the case of the A 
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FIG. 15: Chiral fits to the A ++ '~ and A+'° mass us- 
ing Eq. (28) taking 0,3=3.9 = 0.0855 fm and 0/3=4.05 = 
0.0667 fm determined from Filled triangles show re- 
sults on a L a — 2.1 fm. The result on the 2.7 fm volume 
at /3 = 3.9 is shown with the filled circle. 
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FIG. 16: Chiral fits to the A ++ '" and A+'° mass using 
Eq. (28) with the lattice spacing determined from the 
nucleon mass. The physical point shown by the asterisk 
is not included in the fits. 



The leading one-loop HBxPT result in the case of the A mass has the same form as that for the nucleon mass and 
is given by 



Acim^ - 



25 3H A o 

81 327T/2 - 



(28) 



where m A , is the A mass at the chiral limit and c\ now denotes the coefficient of the m^-term for the A mass. For 
the A axial coupling, Ha, we use the SU(6) relation Ha — (9/5)gA and therefore the one- loop contribution takes 
the same numerical value as in the nucleon case. We also consider the 8— scheme to order C(p 4 /A) which yields an 
expression that is similar to the nucleon case: 
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FIG. 17: Continuum limit of the A + '° mass using the 
lowest order HBxPT to interpolate except for the heavi- 
est pion mass where we use linear interpolation. 
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FIG. 18: Simultaneous chiral fits to the A ++ '~ and A+'° 
mass after extrapolating to the continuum limit at fixed 
rom w excluding the heaviest pion mass. The asterisk 
denotes the physical A mass using the ro = 0.444(3) fm 
extracted in the pion sector and it is not included in the 
fits. The cross denotes the physical point when we use 
the value of ro extracted from the nucleon sector with 
only statistical errors. 



The fits using the 0(p 3 ) HBxPT result at (3 = 3.9 and (3 = 4.05 are shown in Fig. 15 for the A++' - and A+<° 
masses using the lattice spacings determined from / T . It is useful to chirally extrapolate the A mass to see how 
close current results are to A(1232) taking the lattice spacings as determined from the nucleon mass. The fits in this 
case are shown in Fig. 16 again using the 0(p 3 ) HBxPT result. Agreement with the experimental value of the A 
mass is better when one uses the lattice spacing determined from the nucleon mass. This indicates that for baryonic 
observables is favorable to use the lattice spacing determined from the nucleon mass. We give the values of the 
parameters that we extract in Table VI. 

The continuum extrapolation is carried out as in the nucleon case. We show in Fig. 17 the results for the three 
different lattice spacings. As in the case of the nucleon, the continuum limit found by averaging results from (3 — 3.9 
and P = 4.05 is consistent with results at (3 = 3.8. Furthermore, we find that in the continuum limit A ++, ~ and 
A +, ° are degenerate within errors, a result that corroborates absence of isospin breaking. We therefore perform 
simultaneous fits to both A ++,_ and A + '° mass using our continuum limit results at the five smallest pion reference 
masses. These fits using leading chiral perturbation theory and the 8— scheme are shown in Fig. 18. The physical 
point is again not included in the fits. We find a A mass at the physical point that is very close to experiment. Again, 
this agreement improves when we use the value of ro fixed from the nucleon mass. The values that we find for 
and ci from these simultaneous fits using the continuum results are in good agreement with the values extracted for 
A ++, ~ and A +, ° at finite lattice spacing. This points to small cut-off effects and to isospin breaking effects that are 
smaller than statistical errors. 
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VI. CONCLUSIONS 

Using dynamical twisted mass fermions we obtain accurate results on the nucleon mass for pion masses in the range 
of 300-500 MeV. The quality of these results allows a chiral extrapolation using heavy baryon chiral perturbation theory 
to 0(p 3 ). There is agreement among all approaches for this lowest order result. Performing a simultaneous fit to our 
results at the two finer lattice spacings we find a value of the nucleon mass of 0.963 ± 0.012(stat.) ± 0.008(syst.) GeV 
where the systematic error is the difference between the mean values obtained at 0(p 3 ) and 0(p 4 ) HBxPT. Comparing 
results at our three (5- values and at the continuum limit we confirm that cut-off effects are small. Given that this 
leading one-loop result in HBxPT yields good fits to our lattice data we use it to extract the lattice spacing from 
the nucleon mass at the physical point. We find £1,3=3.9 = 0.0889 ± 0.0012(stat.) ± 0.0014(sysi.) fm and 0,3=4.05 = 
0.0691 ± 0.0010(stat.) ± 0.0010(syst.) fm. Again the systematic errors are estimated by comparing the value obtained 
at lowest order to the results obtained using the 0(p 4 ) in HBxPT. Within this estimated uncertainty of the chiral 
extrapolation, the value we find for the lattice spacings 0/3=3.9 and 0/3=4.05 is consistent with the value determined 
from f^. A combined analysis of data in the pion and nucleon sector is a promising option that will be considered 
in the future. We use continuum extrapolated results to determine also the value of r using the nucleon mass at 
the physical point. We find a value of r = 0.473 ± 0.09(stat.) ± 0.016(syst.) fm which is, within errors, consistent 
with the value determined from /„- . The confirmation that isospin breaking in the A is consistent with zero is a very 
important conclusion of this work. This is demonstrated by evaluating the mass splitting in the A isospin multiplcts 
for three lattice spacings on two volumes. Consequently the mass of the A ++, ~ and A + '° obtained in the continuum 
limit are the same within statistical uncertainties. 

The reliable determination of the lattice spacing from the nucleon mass as well as the fact that isospin breaking is 
consistent with zero for these lattices paves the way for further applications of twisted mass fermions in the baryon 
sector. 
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